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Using renormalization-group techniques we analyze equilibrium properties of a large gated quan- 
tum dot coupled via a long and narrow channel to a reservoir of electrons. Treating the electrons 
in the channel as one-dimensional and interacting, we demonstrate that for nearly-open dot and 
not very strong spin fluctuations the ground-state properties as a function of the gate voltage are 
non-analytic at the points of half-integer average dot population. Specifically, the exact result of 
K. A. Matveev, Phys. Rev. B 51, 1743 (1995), that the dot capacitance shows periodic loga- 
rithmic singularities is rederived as a special case corresponding to non-interacting electrons. We 
demonstrate that this conclusion also holds in the presence of SU(2) spin symmetry, and argue 
that logarithmic singularities persist as long as the Coulomb blockade is destroyed which will be 
the case for sufficiently large tunneling and not very strong interparticle repulsions. We show that 
interparticle repulsions aid the Coulomb blockade to survive disordering effect of zero-point mo- 
tion provided the tunneling is sufficiently weak. Upon increase of tunneling, the Coulomb blockade 
disappears through a Kosterlitz-Thouless phase transition characterized by a dot population jump 
whose magnitude is only determined by interparticle repulsions in the channel. Similar conclusions 
are applicable to systems without SU (2) spin symmetry except that the logarithmic singularity of 
the capacitance is replaced by a power-law non-analyticity whose exponent characterizes the degree 
of spin fluctuations. 
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I. INTRODUCTION 

Quantum dots are mesoscopic objects whose proper- 
ties, like those of ordinary atoms, can vary appreciably 
upon addition or removal of a single electron 1 -. This fact 
makes them an active playground of studying the many- 
body physics under the conditions of confinement. 

A simple picture of a dot pertinent to understanding its 
equilibrium properties is that of a container of electrons 
connected to a lead via a single constriction and placed 
nearby a gate electrode. The role of the latter is to tune 
the chemical potential of the dot to host a desirable num- 
ber of electrons which is generally not integer. On the 
other hand, the electrons are discrete, and can only hop 
onto or off the dot one at a time. Therefore they enter 
or leave the dot only at special values of the gate voltage 
corresponding to half-integer dot population. This effect 
well identified by the early 1990's is called the Coulomb 
blockaded. One of its hallmarks is the dependence of the 
population of the dot on the gate voltage consisting of 
plateaus corresponding to an integer number of electrons 
on the dot - the famed Coulomb staircase. 

Quantum-mechanical tunneling of charge between the 
connecting lead and the dot has been argued to smear the 
classical Coulomb staircase even at zero temperature^. 
This idea is supported by a connection of the quantum 
dot problem to a Kondo problem 4 -. In the limit of almost 
zero tunneling the population dynamics on the dot can be 
approximated by that of a two-level system, with the lev- 
els corresponding to population states of the dot differing 
by unityA Changing the population on the dot is then 
akin to flipping an impurity spin-1/2 while the deviation 



of the average dot population from the half-integer value 
is analogous to a magnetic field applied to the Kondo 
spin. For a single-channel Kondo problem pertinent to 
the spinless version of the quantum dot problem, the spin 
susceptibility is finite at zero applied field. This implies 
smearing of the population jumps with variation of the 
gate voltage. 

Complementary treatments of the strong-tunneling 
regime^ led to a conjecture that the rounded Coulomb 
staircase of the weak-tunneling limit smoothly evolves as 
tunneling is increased into a strictly linear function of the 
gate voltage. 

A spinless version of the Coulomb blockade problem 
has been recently revisited^ where additionally it was as- 
sumed that the dot is connected to a one-dimensional 
reservoir so that interparticle interactions have to be 
taken into account. It has been demonstrated that as 
the strength of the tunneling is tuned, changing the dot 
from open to closed, a modified Coulomb staircase reap- 
pears under certain circumstances. This was arrived at 
by mapping the zero-temperature dynamics of interact- 
ing electrons to the classical statistical mechanics of a 
one-dimensional Ising model with ferromagnetic interac- 
tions decaying as the inverse-square of distance. This led 
to the conclusion that all experimental systems must be 
either of the Kondo subclass (to which the connection 
to the Kondo problem 4 applies) or of a novel tricritical 
type. 

Here we reconsider the problem of the Coulomb- 
blockaded quantum dot and take into account both the 
electron spin and electron-electron interactions in the 
channel connecting the dot with the electron reservoir. 
This is relevant because the strongest argument in fa- 
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vor of the Kondo analogy has been given for the non- 
interacting version of this problem: Matveev& found ex- 
actly that in the high-transparency limit, the capacitance 
of the dot has a logarithmic singularity whenever the av- 
erage population of the dot is half-integer. The weak- 
tunneling analysis of this problem^ which used the so- 
lution of the two-channel Kondo problem also found a 
logarithmic divergence which led to the contention that 
this singularity of the capacitance persists for arbitrary 
tunneling. These predictions 6 - found partial experimental 
support^ thus making the issue practically relevant. 

The present work investigates the problem in a phe- 
nomenological fashion by using renormalization-group 
techniques which have an advantage of generality, thus 
complementing earlier calculations. Specifically, we show 
that the logarithmic singularity of the capacitance^ is a 
consequence of the assumption that the electrons in the 
constriction are non-interacting. We also point out that 
the same conclusion holds in the presence of SU (2) spin 
symmetry as long as the Coulomb blockade is destroyed. 
Moreover, renormalization-group methods allow us to see 
clearly possible limitations of existing approximations, 
thus motivating future work. 

The organization of this paper is as follows. In Sec- 
tion II we set up the problem in the form of a zero- 
dimensional quantum field theory describing the dynam- 
ics of the charge and spin of the dot. For a nearly-open 
dot (Section III), the population of the dot tracks the 
gate voltage and only the spin dynamics described by a 
sine-Gordon like theory is relevant. In Section IIIA this 
field theory is studied perturbatively in small reflection 
amplitude. This exercise demonstrates that physically 
interesting cases are beyond the perturbative approach. 
Section IIIB contains a heuristic argument which allows 
us to guess the answers for situations when the pertur- 
bation theory fails. These conclusions are placed on a 
sounder footing in Section IIIC where renormalization- 
group analysis is employed. The novel feature here is 
a systematic calculation of the renormalization of the 
ground-state energy which allows, for example, to exactly 
reproduce the logarithmic divergence of the capacitance 6 . 
In Section HID we draw a parallel between these results 
and those found earlier for the standard sine-Gordon the- 
ory. 

In Section IV we study the opposite limit of a nearly- 
closed dot. For half- integer average dot population the 
dynamics of electrons hopping onto and off the dot is 
approximated by that of a two-state system. The zero- 
temperature properties of the latter can be inferred from 
the classical statistical mechanics of a two-component 
model resembling the one-dimensional Ising model with 
interactions decaying as the inverse-square of distance. 
An "entropy" argument which in Section IVA is applied 
to this model provides us with the range of stability of 
the regime of zero tunneling. These conclusions are re- 
fined in Section IVB via a renormalization-group calcu- 
lation perturbative in tunneling amplitude. This calcu- 
lation indicates that whenever the tunneling is relevant, 



the non-analytic behavior of the ground-state properties 
of the dot has to be captured by the results found for a 
nearly-open dot. In Section IVC this statement is made 
more quantitative via a heuristic argument which shows 
how the weak-tunneling results away from half-integer 
dot population can be inferred from those of Sections 
III and IVB. We conclude (Section V) by summarizing 
our results and emphasizing possible limitation of the 
two-state approximation of the hopping dynamics for a 
nearly-closed dot. 



II. STATEMENT OF THE PROBLEM 

Following Matveev 6 - we consider a large quantum dot 
of capacitance Co whose charging energy Ec — e 2 /2Cq is 
significantly larger than the distance between the energy 
levels on the dot. The dot is placed near a gate and con- 
nected to a reservoir of spin- 1/2 electrons through a chan- 
nel which we assume to be so narrow and long that the 
electrons in the channel are essentially one-dimensional 
and their interactions have to be taken into consideration. 
Although there are no known experimental realizations of 
this model, appropriate experiments might become pos- 
sible in the future. Additionally, this represents a gen- 
eralization of the model due to Matveev 6 - in a way that 
adds interactions violating SU(2) spin symmetry because 
it is interesting from theoretical point of view. The cor- 
responding long-wavelength low-energy Euclidian action 
has the form 



1 



1 = -/ ^d X dr Pu U^r 

iy—c,s 



dr7(u T (0,r),« i (0,r)) 



4<£)" 



(1) 



The first harmonic term describes interacting one- 
dimensional electrons inside the channel. Here the dis- 
placements of the electrons of each spin population with 
respect to their classical positions are represented by the 
fields u^(x, t) and r) where arrows indicate the spin 
direction, x is the one-dimensional coordinate along the 
channel and r is the imaginary time. Written in terms 
of the center of mass ("charge") and reduced relative dis- 
placement ("spin") fields 



(2) 



correspondingly, the first term of ([I]) is diagonal. The 
parameters p v and c„ are the macroscopic mass densities 
and sound velocities in the charge and spin sectors^. 
The harmonic (Luttinger) liquid exponents 
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characterize interparticle interactions in the charge (v = 
c) and spin [y = s) channels^. They are formed out of 
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the parameters of the harmonic part of the action ([I]) 
and the one-dimensional electron density n within each 
spin population. Although non-interacting electrons cor- 
respond to g c , s — 2, in what follows we will consider more 
general values of the exponents. 

We also note that either spin population creates a pe- 
riodic potential for its opposite spin counterpart which 
may cause spin ordering. For weak short-range inter- 
action between the electrons of opposite spin this can 
be modeled by adding to HJ a term proportional to 
J dxdr cos 47m it,,.. We however assume that such interac- 
tions are irrelevant in the renormalization-group sensed, 
and thus the spin ordering does not take place. This 
would be the case, for example, in the presence of SU (2) 
spin symmetry (the continuum Hubbard model with 
repulsive interactions can serve as an example) where 
g s = 2. More generally we will require that the degree of 
quantum fluctuations in the spin channel is sufficiently 
large^, specifically g s > 2. 

The second term of the action ([1]) is due to the tunnel- 
ing and Coulomb repulsion on the dot and the arguments 
of the potential V are evaluated at x = which can be 
taken to be the narrowmost place of the channel^. 

Since the interactions in ([1]) are restricted to a point, 
the harmonic degrees of freedom away from x = can 
be integrated out. This is equivalent to minimizing the 
action {TJ with fixed values of the fields u CjS at the origin 
followed by the substitution of the solution back into ([I]) 
and evaluation of the outcome. The conditions 5A/5u c , s 
= reduce to the Laplace equations 

whose solutions satisfying the boundary conditions 
w CjS (0, r) = u c , s (t) have the form 



duj 
2^ 



u CjS (ui) exp(iarr - |w||a;|/c c , s ) (5) 



where we introduced the Fourier transforms u c . s (uj) = 
J dru CiS (r) exp(— iuir). Substituting © in ([J), perform- 
ing the integration, and introducing two new fields 9 C 
and 9 s /2 corresponding to the fluctuating electron num- 
ber and spin on the dot, respectively 

9 C = n(u-\ + u\) = 2nu c , 9 S = n(u-\ — ttjj = 2nu s (6) 

we find the effective action 

' Lj M|0 w (u;)| 3 + J drV{9 c ,9 s ) 
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(7) 

where Ao is the high-frequency cutoff. Here without the 
loss of generality we restrict ourselves to the sharp cutoff 
model with all the \uj\ > Ao modes eliminated. 

For the problem at hand the potential V entering the 
actions ([1]) and (7]) is a sum of two contributions: 

- Nf + U(u hUl ) = E C (9 C - N) 2 + U(6 C , 6.) 

(8) 



V = Er 



where the first term is the charging energy with dimen- 
sionless parameter N proportional to the gate voltage 6 . 
The second, tunneling term U{w[,U\) is periodic in both 
arguments with period n~ 1 which reflects the fact that 
as either of the spin populations moves along the channel 
it experiences a periodic potential with the period set by 
interparticle spacing. Equivalently, it is a manifestation 
of the discreteness of the electrons - the dot "prefers" to 
host an integer number of electrons thus implying that 
the function J7(iif, ujj has minima at integer values of 
nuf t i. The classical limit of "closed" dot is recovered 
by increasing the amplitude of U to infinity while the 
regime of perfect transmission is reached by decreasing 
the amplitude of U to zero. For example, for the reflec- 
tion amplitude satisfying the condition \r\ 1, the form 



D\r\ 

U = (cos 2irnui + cos 27mwj) 

7T 

27A ft|r| 
= costo,- cos Try, 



(9) 



was derived earlier—. Here D is the high-energy cutoff 
(bandwidth) used by Matveev 6 - and 7 = with C w 
0.5772 being Euler's constant. Shortly it will become 
clear why the energy and frequency cutoffs in |(9]) are 
related by D = frjAo. 

We note that without the charging term, Eq = 0, 
Eqs.® and describe the dynamics of spin-1/2 elec- 
trons in the presence of a localized inhomogeneityii. Ad- 
ditionally, for the special case g c = g s these equations 
capture the physics of resonant tunneling of interacting 
spinless electrons onto and off a dot formed by two tun- 
neling barriers in the presence of a gate which couples to 
the electrons between the barriers^i. 



III. NEARLY-OPEN DOT 

Combining Eq.((7]), ([8]), and ([9]) we arrive at the effec- 
tive action describing the physics of nearly-open dot 
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(10) 

For almost perfect transmission \r\ <C 1 the action (jTOJ) is 
minimized for 9 C = N , and the electron number field 9 C 
can be eliminated. This is done by writing 9 C = N + 59 c 
where S9 C is a fluctuation and averaging over the unper- 
turbed, \r\ — 0, ground state, i. e. over the Gaussian 
part of (fTTj|) involving the fluctuating number field 9 C : 



duo . irh . _ , . . , , . o 
—(—\ U \+E c )\6 c (u)\ 3 11 

Z7T g c 
|w|<A 



Aoc 
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This can be used to calculate the mean-square fluctuation 
of the number of electrons on the dot 
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where <>o c denotes an expectation value computed with 
Aoc, Eq. (fTTj) . Similar to the spinless case 5 -, the fluctua- 
tion lfl~2j) is finite due to the presence of the Coulomb re- 
pulsion on the dot. We note however that for the special 
case of free electrons, g c — 2, Eq . (|12[) agrees with pre- 
vious conclusion 6 only with logarithmic accuracy. When 
the frequency cutoff A vanishes or the charging energy 
Ec diverges, the fluctuation lfl2|) vanishes. Additionally, 
when Ec — 0, the fluctuation (fT2]) diverges. These con- 
clusions are in accordance with the physical expectation. 
The result lfl2"Tl allows us to find the expectation value 



< cos7T0 c >0c= (1 



Ec9c' 



^ /4 cos7rN (13) 



With this in mind the action (jTUJ) simplifies to a reduced 
form involving only the spin field 6 S 



A s f dm 

=77 / 

h J 2ng t 

M<Ao 



\uj\\8 s (uj)\ 2 -v A / dr costtBs (14) 



where the dimensionless parameter 

,o = ^(l + ^r- /4 co S7 ^ (15) 
7T E c g c 

accumulates salient information about the problem at 
hand, and the parameter TV has a meaning of the av- 
erage population of the dot. 



A. Perturbation theory 

We start by computing a correction to the ground-state 
energy of the dot SE perturbative in the interacting part 
of the action (Tl4l) 



dr cos 7T#, 



(16) 



A quick way to arrive at corresponding expression is to 
notice that the action (fl4|) can be formally recognized as 
a Hamiltonian of a one-dimensional classical statistical 
mechanics problem with the imaginary time coordinate 
t corresponding to fictitious space directio n 12 ' 13 . Then 
the zero-point motion is analogous to thermal fluctua- 
tions, and the result can be inferred from the formulas of 
thermodynamic perturbation theory^. Restricting our- 
selves to second order we have 

L5E =< Ai n t >Q S — — < (Ai nt — < Ai n t >0s) 2 >0s 

(17) 

Here L — h/T is the size of the fictitious classical system 
in the imaginary time direction and T is the temperature; 



the zero-temperature limit is equivalent to the L — > oo 
limit. The notation <>q s indicates an expectation value 
computed with the harmonic, v = 0, part of A s , Eq.fl"4"|): 



A 0s 

h 



l<Ao 



dw 
2ng s 



(18) 



In order to compute the first-order correction, 8Ei =< 
Am >0s / L, we need to know the mean-square fluctua- 
tion of the spin field 9 S which can be found to be 



^ nl ^ _ Us 



|w|<A 

This implies that 

LSEk =< Ant >0s= -voHAoL) 1 -^/ 4 , (20) 

i.e. the first-order correction to the ground-state energy 
SEi oc L~ 9a / 4 vanishes in the zero-temperature limit. 
Similar to the free-fermion case£, g c , s — 2, this is a con- 
sequence of divergent spin fluctuations lfT9|) . 

In order to calculate the second-order correction SE 2 = 
- < (A int - < Ant >o s ) 2 >o s /2LH, we will need the 
correlation function 



X (r) =< (6 s (t) - 9 S (0)) 2 > 0s 
g s f du 1 — cos lot 

7 J 2^ M 

\w\<A 



r Cin(A |r|) (21) 



where C'm(t) is a cosine integrate. As t — > the func- 
tion Cin(i) vanishes as £ 2 /4 while as t — * oo, it diverges 
logarithmically, Cin(i) — » ln(7t). The latter asymptotics 
also implies the large-argument logarithmic divergence of 
the correlation function l[2"Tj) . x(r) = (g s /ir 2 ) ln(7Ao|r|). 
For noninteracting electrons, g s = 2, this can be com- 
pared with the result implicit in Eq.(34) of Matveev^: 
x(t) = (2/-7T 2 ) \n(D\r\/h). Therefore the energy and the 
frequency cutoffs, D and Ao, respectively, are related by 
D = hjAo, which was already employed in Eq.®. 

The expression for the second-order correction 8E2 can 
be simplified to 



SE 2 = - 



2h 



00 

J drexp(- 



where the function 



B(g s ), 
(22) 



OO 

B{g.) = /dtexp(-|Cin(t)) (23) 



is only defined for g s > 2. The evaluation of the integral 
in Eq. (|23|) employing the logarithmic approximation of 
the cosine integral, Cin(f < 7 _1 ) = and Cin(t > 7 _1 ) = 
ln(7*), gives 



B(g s ) 



7(g s - 2) 



(24) 



5 



Although this estimate becomes asymptotically exact as 
g s — > 2 from above, as g s increases, its accuracy degrades. 
Specifically, as g s — > oo, the integral lf23|) is dominated by 
the region of small f where the logarithmic approximation 
is not adequate and we find instead B(g x ) « (27r/g , s ) 1 ' 2 . 
However exceedingly large values of the spin correlation 
exponent <? s are non-physical and will not be considered. 

In order to evaluate a low-temperature contribution 
into 5E 2 we need to replace the upper integration limit 
in Eq.([23| with A L/2 > 1. As a result we find 



SE 2 = 



VqHA 



B(g s ) 1 



B(g s 



2( 7 A L/2) 1 -^/ 2 
7 ( 5s - 2)B(g s ) 
2( 7 A L/2) 1 -^/ 2 



(25) 



where in the last step we used the approximation (|24|l . 
For g s > 2 we may take the limit of zero temperature, 
L — > oo, which takes us back to Eq. (|22|) . while in the 
presence of SU(2) spin symmetry, g s — 2, we find with 
logarithmic accuracy 



5E? 



2 7 



In(AoL) 



(26) 



In the zero-temperature limit, L — > oo, this result di- 
verges which means that a non-perturbative calculation 
is required. The logarithmic divergence f26|) is analogous 
to its free-fermion counterpart g CtS — 2 found earlier—. 

Additionally, the special place of the limit of non- 
interacting electrons can be seen by taking the contin- 
uum limit, Ao — » oo, in Eq. f26|) . This only leaves 
the logarithmic dependence on the microscopic cutoff, 
SE 2 = -(47|r| 2 £;c/7r 3 )ln(Aoi), thus implying that the 
prefactor is solely determined by macroscopic quantities. 



Beyond perturbation theory: 
argument 



heuristic 



In order to go beyond the perturbation theory, first 
we need to clarify its range of applicability. A naive ex- 
pectation that the perturbative result lf22]) is reliable for 
g s > 2 as long as \v$\ <C 1 is misleading. To find the 
correct condition we adopt an argument due to J. La- 
jzerowicz originally given in the classical context of the 
roughening phase transition^: 

The cumulant expansion (fT7|l is performed in powers 
of Aint/h which implies that the perturbation theory is 
applicable as long as | < A int >o s <h Inspection of 
Eq, (|20p shows that this will be the case only for g s > 4 
and | | <C 1. One can also arrive at the same conclusion 
by applying the standard argument which compares 
the magnitudes of the first and second-order terms of the 
expansion (|17p . 

The physics in the non-perturbative regime 2 < g s < 4 
can be understood heuristically by noticing that the per- 
turbative treatment assumes that spin fluctuations grow 



with scale without bound (see Eq. lfT9|) ) which is indeed 
the case for g s > 4. This assumption, however, breaks 
down for g s < 4 because spin fluctuations are suppressed 
at scales exceeding a characteristic value L c which has a 
meaning of the correlation length of the fictitious clas- 
sical statistical mechanics problem corresponding to the 
action <fl4|l . This scale marks the verge of applicability 
of the perturbation theory, | < Ai nt /h\ ~ 1, which with 
the help of Eq. (|20|) leads to the estimate 



An 



-4/(4- s .) 



(27) 



In the limit g s — ► 4 from below the correlation length 
(|27ll diverges which is in correspondence with the well- 
known factA£ that the model lfl4|) has a phase transi- 
tion at g s — 4. We also note that for the case of non- 
interacting electrons, g c . s — 2, the energy scale h/L c cor- 
responding to lf27j) has been identified earlier—; a related 
renormalization-type argument is given in Section IIIC. 

At the scales exceeding L c the spin fluctuations are ef- 
fectively frozen and no longer contribute into the ground- 
state energy. Thus the latter can be estimated by sub- 
stituting L ~ L c in Eqs. (|25|) and (|26| . As a result for 
2 < g s < 4 the correction to the ground-state energy SE 
is estimated as 



5E~ - 



B(g s ) 1 



\vq 



2( S .-2)/(4-fl,) 



while for g s = 2 we find 



SE 



Vgfi A 
2 7 



In 



(28) 



(29) 



The numerical uncertainty present in the definition of 
the characteristic scale L c , Eq. lj27|l . translates into an 
uncertainty of the term of Eq. l(28|) non-analytic in the 
parameter vq: here vq should be viewed as vq times a 
constant of order unity. However the analytic Vq term 
is given exactly. At the same time Eq. (|29|) continues 
to possess logarithmic accuracy as the uncertainty of L c 
only affects the argument of the logarithm. 

Since the parameter vq, Eq. (fl~5|) . vanishes at half- 
integer dot population N, the ground-state energy given 
by Eqs. ((2Sj) and (|29|) and related equilibrium prop- 
erties exhibit non-analytic behavior at half-integer N. 
For example for non-interacting electrons, g c , s — 2 (or 
more generally in the presence of SU (2) spin symmetry, 
g s = 2), the capacitance of the dot C oc d 2 E/dN 2 di- 
verges logarithmically at half-integer dot population, the 
result derived by Matveev through an exact calculation^. 

More generally, 2 < g s < 4, the singular part of the 
energy behaves as 8E slng oc \SN\ 4 ^ 4 ~ 9 ^ (see Eq.pS])) 
where SN is a small deviation from a half-integer dot 
population. Although the corresponding contribution 
into the capacitance C S i ng oc \SN\ 2 ^ 9 "~ 2 ^^ 4 ~ 9 "^ vanishes, 
the third-order derivative of the energy, d 3 E sing /dN 3 oc 
\§ N \(3g.-8)/(4-g.) diverges for g s < g/ 3> Generally, the 
spin correlation exponent g s parameterizes the order of 
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the phase transition at half-integer population of the dot 
(i.e. the order of the first divergent derivative of the 
ground-state energy): it is of the second-order at g s = 2 
and of the Zth-order for g s < 4(1 - I) /I and I > 2. As 
g s increases from 2 to 4, spin fluctuations increase, the 
non-analyticity of the ground-state energy manifests it- 
self in the divergence of progressively higher-order deriva- 
tive and the order of the phase transition increases ap- 
proaching infinity at g s = 4. The ultimate origin of these 
singularities is the discreteness of the electron spin. For 
g s > 4 spin fluctuations become unbound and the cor- 
rection to the ground-state energy becomes regular and 
proportional to -(cos2?riV + 1) (see Eqs.([22]) and [fl5j) ) 
still exhibiting remnants of the discreteness of the elec- 
trons. 



C. Beyond perturbation theory: 
renormalization-group analysis 



The heuristic results (|28j) and (|29|l can be placed 
on a firmer ground through a perturbative, »o < 1, 
renormalization-group analysis of the effective action 
(jT4j) . It is well-known that renormalization procedure 
involves repeated elimination of microscopic degrees of 
freedom^. This process generates terms which renor- 
malize the parameters of the initial action thus bringing 
about observable physics, and contributions which are ir- 
relevant in the renormalization-group sense and thus can 
be dropped. For the model (fl4l) the associated terms are 
discussed in the literature^. Additionally, elimination 
of microscopic degrees of freedom can generate exten- 
sive quantities which (in the present context) represent 
contributions into the ground-state energy. As far as we 
know, these contributions which will be our primary con- 
cern, have not been discussed which gives us an excuse 
to summarize the salient steps of the whole calculation. 

In the sequel we adopt the approach of Nozieres^ 
who investigated the classical problem of the roughen- 
ing phase transition. As usual, the field S is presented 
as a sum of "fast", 9P, and "slow", 9 S , modes: 



S ( T )=0(/)( T )+0W( T ) 



(30) 



Below we will use the sharp cutoff procedure^ appropri- 
ate for the model lfl4|) . i. e. choose 



(31) 



0P(w) = 6.(w), A (l-e)<H<A 0> 



where e <C 1, and 



M < Ao(l-e) 



(32) 



After the fast degrees of freedom are eliminated, the ef- 
fective action (f!4|) turns into a functional of slow modes 

a s [oP 

expansio 



r)j, which to second order in vq is given by the 

16.17. 



M<A (l-e) 



2ng s 



+ <^, t [fiW(r)+^(r)]>/ 
- ^<{Aint\6 s a \r) + 9{f\r)] 
- <A mt [9i s Hr)+0i f Hr)}> f ) 2 > f 



(33) 



Here <>/ denotes averaging over fast degrees of freedom 
with the Gaussian probability distribution determined by 
the action lfT8|) except that the frequency range, accord- 
ing to Eq.([3T]). is A (l - e) < \uj\ < A . 
To first order we find 

Aj = -v (l - y < (^P? >f)Mh J drcosnOP 
= -v (l - je)A ft J drcosnOP 



vAh / dr cos n9i s \ 



(34) 



where we employed the fact that the mean-square fluc- 
tuation of the fast modes, < (0^) 2 >f, is given by 
the integral lfT9|) except that the integration range is 
Ao(l - e) < \uj\ < A , Eq.((3l|). The last representa- 
tion in (|34|) combined with the first term of <f33|l form an 
action having the same functional form as the original 
action A s , Eq. (fl4|) except for smaller cutoff, 



A = A (l-e), (35) 
and a slightly different amplitude of the cosine term, 

i - 9se/4 



v = v - 
To second order we obtain 



1 -e 



(36) 



A-2 



drdr' sin7r6»i s) sin7r6»i s ' 



x <0^9^'> f , (37) 



where the prime means that the argument of corre- 
sponding function is r'. The correlation function < 
ei f) ei f) ' >/ is given by 



< op op 



>f 



du expiw(r — r') 



9s 

2rr J 27: 

Ao(l-d;)<l"l<Ao 

Jk 

2rr 2 



-ecos(A (r - r')), 



(38) 

which allows us to rewrite Eq. f37|) in a more explicit form: 

v 2 Alhg s 



Ao 



e / drdr' cos(A (r - r')) 



x (cos7r(^ s ) + OP ') - cos7r(6l( s ) - OP '))(39) 

A further truncation is necessary to see the implica- 
tions of Eq.j39|- First, the cos n(oP +oP ') term can be 
dropped since for r = r' it becomes the second harmonic 
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of the main cos ir9 s oscillation, and thus is irrelevant in 
the renormalization-group sense. 

Gradient expansion of the other term 

cosTr^W-flW ') « l-f.(0[')-6W >f 



« l-y((r-r')^ s )/dr) 2 (40) 
is not appropriate for large magnitude of the difference 



t— t' where the two values 9s and Of ' are hardly corre- 
lated: the right-hand side of (|40)l needs a factor specifying 
the extent to which these values are correlated. The same 
situation is encountered in the context of the usual sine- 
Gordon theory where Knops and den Oudeni^ demon- 
strated that continuing the expansion of the cosine to all 
orders generates the necessary factor. Below we adopt 
an equivalent physical argument due to NozieresA^, and 
go beyond the existing treatments of the action (|14p : 
The field 9^ in (|40)l should be considered as a sum 

(s) 

of the equilibrium fluctuation 9 S and a small correction 
C (which might be a response to a small low-frequency 
perturbation) . Thus instead of (|40l) , first we average over 
equilibrium fluctuations which is followed by the expan- 
sion in C- This generates the expansion replacing Eq. ij40|l : 

costt(0W -9 ( s s) ') 

C') 2 )<COS^-^ ')>„, 



5« / 



1 



((T-T , )d9^/dl 



x exp(- 



2 

tt 2 x(t 



(41) 



We note that the same correlation exponential appears in 
the second-order result, Eq. ([22|) . >From technical view- 
point, the presence of this factor eliminates all the po- 
tential difficulties usually encountered when sharp cutoff 
procedures are use d 10 ' 12 . 

Insertion of the expansion l[4"T]) back into |39|) makes 
it clear that a small positive contribution of the 
/ dr{d9 s /dr) 2 form is generated under renormalization. 
For the classical problem of the roughening phase transi- 
tion such contributions are important parts of the physics 
as they renormalize the interface stiffnessAS'. However in 
the present context they are irrelevant and thus can be 
dropped^. For the model lfl4|) the role of the stiffness is 
played by gj 1 which does not renormalize to any order—: 
renormalization procedure cannot generate non-analytic 
contributions proportional to \uj\. 

The constant term in (|4l"Tl generates a small negative 
correction to the ground-state energy 



where the function 



v 2 A h 



H(g a )e, 



X 

H(g s ) = y / dt cost exp(~^Cm (*)) 



(42) 



(43) 



is defined unambiguously for g s > 2. Using the definition 
of the cosine integral l|2"Tj) this formula can be rewritten 
as 



H{gs) = |/(l-*^»)exp(-|Cin( W 



t i-gs/2 d exp (_^ Cin ( i ))^ (44) 



which implies 

H{2) = lim (texp(-Cin(t))) = 7 _1 w 0.5615 (45) 

t — >oc 

This demonstrates the special role played by g s = 2 - 
only in the presence of SU(2) spin symmetry is the cor- 
responding value of H(g s ) controlled by the t — > oo limit. 

For g s > 2 we integrate Eq. (|44j) by parts which es- 
tablishes a relationship with the function B{g s ), Eq. ([23l) . 
appearing in the second-order expression (|22| : 



H(g s ) = (| 



1)-B( 5 a) 



(46) 



Renormalization-group transformation which led to 
Eqs. (|35)l . (J36j) , and (|42|) can be carried over repeatedly 
by reducing the frequency cutoff by successive infinitesi- 
mally small amounts e = dl — > 0. Since each step of the 
transformation incorporates the effect of previous renor- 
malizations, the outcome is customary to represent in the 
form of differential flow equations: 



dk 
~dl 



-A 



dv 
~dl 



I 1 - -r) v 



dE 

~dJ 



v 2 Ah 



H(g s ) 



(47) 



(48) 



(49) 



These equations have to be solved subject to the "initial" 
conditions, A(l = 0) = A , v(l = 0) = v ; the macro- 
scopic behavior is recovered in the A(l — ► oo) — > limit. 
Eq. (|47| whose solution is given by 



A(0 = A e 



(50) 



is a common place in renormalization-group calculations. 
Eq.jHJ and its solution 



v(l) = v e^- g ^ 1 



(51) 



are also well-known^ from renormalization-group treat- 
ments of the model lfT4|) : 

For g s > 4 the parameter v(l) flows to zero as I — » oo 
thus implying that the last term of lfl4|) is irrelevant in 
renormalization-group sense. This is the region of ap- 
plicability of perturbative approach with unbound spin 
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fluctuations on the dot (see Eq. (fT9|) ). For g s < 4 the 
parameter v(l) grows under renormalization indicating 
that the discreteness of the electron spin is relevant and 
that perturbative treatment fails at low frequencies. As 
a result the fluctuations of the spin of the dot will be 
suppressed at flow parameters exceeding a characteristic 
value l c such as \v(l c )\ ~ 1: 



■ III 



.9* 



1 

N 



(52) 



We note that the corresponding scale in imaginary time 
direction At 1 = An e lc coincides with the heuristic result 



(|27[) - the argument that led to ([27)1 is equivalent to the 
first-order renormalization-group calculation. 

Eq. (|49|) describing renormalization of the ground-state 
energy due to eliminated high-frequency modes is novel; 
its counterpart in the context of the conventional sine- 
Gordon problem is know n 12 ' 16 . In order to compute the 
contribution into the ground-state energy we first substi- 
tute the solutions 1|50|1 and l|51|) in Eq. (|49|) which turns 
it into a form 



dE 



v 2 A h . 



-H(g s )e 



(i-9»/2); 



(53) 



that makes explicit a special role played by SU(2) spin 
symmetry, g s = 2. 

For g s > 4 spin fluctuations of all frequencies con- 
tribute into the ground-state energy thus implying that 
Eq. i|53p should be integrated between zero and infinity. 
Combined with Eq . (|46|) . it then reproduces the pertur- 
bative result (|22|) . 

For 2 < g s < 4 spin fluctuations are suppressed at 
small frequencies and the contribution into the ground- 
state energy can be found by integrating Eq. ([53]) between 
zero and the characteristic flow parameter l c , Eq. (|52|) . As 
a result for 2 < g a < 4 we find 



SE 



VqAqH 



B(g s )(l - \vo\ 2 ^- 2)/{4 -^) (54) 



This is very similar to the heuristic result (|28| . and the 
uncertainty of the second term of (|54|) having its origin 
in the uncertainty built into the definition of the charac- 
teristic flow parameter l c , Eq. ([52"ll . makes Eqs. f28|) and 
(|54| practically indistinguishable. 

Finally, for g s = 2 integrating Eq. ([53"ll between zero 
and l c , Eq. (|52|) . and taking into account the numerical 
value if (2), Eq . lj45|i . reproduces the heuristic result J29|)- 



D. Analogy to the sine-Gordon theory 

The main conclusion of previous analysis is that in 
the free-fermionic/S'[/(2) symmetric limit, g s = 2, the 
ground-state energy (as a function of deviation 5N away 
from half-integer dot population) exhibits logarithmic 
non-analyticity, 5E sing oc 8N 2 \n(l/\8N\); in the inter- 
acting case, 2 < g s < 4, this is replaced by a power law, 



5E S i n g oc \8N\ 4 ^ 4 ~ 9s K This outcome is not unique to 
the zero-dimensional field theory lfl4|) describing the dy- 
namics of nearly-open quantum dot - the conventional 
one-dimensional sin-Gordon field theory shows the same 
effect. The pertinent physical phenomenon here is the 
Peierls instability-^: 

More than fifty years ago Peierls demonstrated that a 
one-dimensional metal is subject to a spontaneous lat- 
tice distortion that converts it into an insulator. The 
distortion causes commensuration of electronic and ionic 
subsystems which leads to the reduction of the electronic 
part of the system energy. For non-interacting electrons 
Peierls found that the energy gain of commensuration as 
a function of small lattice distortion S has non-analytic 
form SE S i n g oc S 2 ln(l/<5). Since other effects favoring 
regular lattice spacing 6 = increase the energy by an 
amount proportional to S 2 , the distortion is inevitable. 

Since the Peierls phenomenon just barely occurs for 
non-interacting electrons, interelectron interactions have 
a drastic effect on the fate of the instability. For spin- 
less electrons the phenomenological theory of the effect 
which takes into account interactions, has been given in 
the past^S. There it was demonstrated that the elec- 
tronic energy gain due to the lattice distortion is given 
by the ground-state energy corresponding to the sine- 
Gordon action 



.4 



A cos lixnu 



(55) 

Here u(x, t) is the electron displacement field, p is the 
mass density, c is the sound velocity, n is the number 
density, and A is proportional to the amplitude of the lat- 
tice distortion 5. The effect of interelectron interactions 
can be parameterized by dimensionless combination 
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irhn 2 
pc 



(56) 



which can be recognized as the "spinless" counterpart of 
the charge correlation exponent g c (J3j) . Renormalization- 
group analysis of the action (|55| resembling that of previ- 
ous subsection then leads to the following conclusions^: 

For non-interacting electrons, g = 1, the form, origi- 
nally derived by Peierls^, SE S i ng oc (5 2 ln(l/(5), is found. 
This is analogous to the "non- interacting" result, Eq. l[29|) . 

In the presence of interactions satisfying g < 2 (and 
excluding g = 1), the singular part of the ground-state 
energy as a function of the lattice distortion S is given 
by 8E sin g oc <5 2 /( 2 ~f). This is analogous to our result 
5E sing oc |<WV| 4/(4 ^ } (see Eqs. (|28f and (|54j)). 

One of the important differences between the actions 
(fl4| and f55|) is that the spin correlation exponent g s 
controlling the degree of zero-point motion in lfl4|) 
does not renormalize. For the sin-Gordon action (|55|) the 
same role is played by the exponent g (|56l) which does 
renormalize. However for g < 2 and small distortion 
(small A) this renormalization is negligible™ That is 
why an analogy can be drawn between the ground-state 
properties of the field theories l|14p and (|55l) . 
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IV. NEARLY-CLOSED DOT 



Although the analysis of Section III makes it clear that 
the logarithmic divergence of the capacitance at half- 
integer dot population^ is a consequence of assumption 
that the electrons in the constriction are non-interacting, 
it cannot explain why the same type of singularity may 
be present in the weak-tunneling regime^ and possibly 
for any tunneling^. To address these issues, we turn our 
attention to the limit of a nearly-closed quantum dot. 

The starting point of the analysis is clarification of the 
nature of the classical ground state described by the ef- 
fective action J7]). This is determined by the minima of 
the potential V(9 C , 9 S ), Eq.®, and has been discussed by 
Kane and Fisher— in a related context of resonant tun- 
neling of spinless electrons confined between two barriers 
in the presence of a gate. Although the potential entering 
Eq. i|10p is not quantitatively correct for a nearly-closed 
dot, it nevertheless retains the symmetries of the prob- 
lem and can be used as an illustration of the argument 
below disregarding the \r\ -C 1 constraint. Visual aid is 
provided by Fig. 1 which parallels Figure 2 of Kane and 
Fishepii. 

Without the charging term, E c = 0, the classical 
ground-state is degenerate with infinitely sharp minima 
located at integer values of the dot population 9 C , which 
reflects the discreteness of the electrons. The charging 
term of iJH) lifts this degeneracy so that for most values 
of the dimensionless gate voltage N there is a unique 
ground state determined by an integer 9 C closest to N 
with energy gaps of the order Ec or larger to the states 
of higher or lower dot population. 

For half-integer values of the dimensionless gate volt- 
age N the ground state is doubly degenerate. Since all 
other integer population states are separated from these 
two by energy gaps, it seems plausible^* that they can 
be eliminated from the discussion of the hopping dynam- 
ics when the dot is slightly open and charge and spin 
fluctuations are included. The validity of this two-state 
approximation will be discussed below but for now we 
assume that it is justified. Then the 8 C coordinates of 
the minima of the effective potential V(0 C , 9 S ) are either 
N + 1/2 or N — 1/2 while their 9 S positions are integers 
with the constraint that for every minimum its 9 C and 
9 S coordinates have the same parity. This means that if 
nearest minima are connected by straight line segments 
(see Fig. 1), they will form an infinite zigzag running on 
average parallel to the 9 S axis. 

Opening the dot means that the barriers separating 
degenerate states of integer population are no longer in- 
finitely high - the effective potential will now have degen- 
erate minima whose 9 C positions are either N + q/2 or 
N — q/2 where q is the distance between the minima of 
V along the 9 C axis^k At the same time the step of the 
zigzag in the 9 S direction which is the symmetry of the 
problem remains unchanged. Including charge and spin 
fluctuations, i. e. allowing finite g c _ s , Eq. ([3]), has the 
same effect of downward renormalization of the distance 



N — 



Figure 1: Positions of the degenerate minima of the potential 
V(6 c ,8 a ) entering Eq.© in the (0 C ,6 S ) plane for half-integer 
dimensionless gate voltage N. For the purpose of illustration, 
nearest minima are connected by dashed line segments. The 
parameter q has a meaning of the distance between the min- 
ima in the 8 C direction. Although q = 1 for a nearly-closed 
dot, zero-point motion renormalizes it to a < q < 1 value. 
In the drawing we assumed (without the loss of generality) 
that N + 1/2 is odd. 



between the minima of V along the 9 C direction. In- 
deed, even for a nearly-closed quantum dot where q = 1 
initially, quantum-mechanical tunneling between nearest 
degenerate ground-states reduces it to a < q < 1 value. 
For q — the effective potential energy evaluated at 
9 C = N is ^-independent. We note that observable value 
of the parameter q corresponds to the jump of the dot 
population as the dimensionless gate voltage N is swept 
across a half-integer value. Then q = means that the 
Coulomb staircase is smeared. 

The elementary act of tunneling between nearest 
ground states consists of a change of the 9 C coordinate 
by q accompanied by simultaneous change of 9 S by unity. 
The zigzag sequence of the ground states implies that 
the tunneling event A9 C — ±q must be followed by an- 
other tunneling event consisting in opposite sign change 
of A9 C = ^fq. At the same time tunneling changes of 9 S 
by unity are not correlated. 

To gain better understanding of the physics in the 
regime of a nearly-closed quantum dot we invoke a 
mapping^ between zero-dimensional quantum field the- 
ory defined by the action ([7]) and that of one-dimensional 
classical statistical mechanics. Similar to the observation 
made for the spinless version of the problem^, Eq. (J7j) can 
be viewed as an effective Hamiltonian of one-dimensional 
classical statistical mechanics problem involving fluctuat- 
ing two-component field 9 C . S with imaginary time r play- 
ing a role of fictitious space direction. The |w| depen- 
dence of the first term of Eq. J7]) implies a 1/r 2 interac- 
tion between the fields 6 CjS separated by distance r. 

If not for the presence of the spin field 9 S , then for 
half-integer N the classical statistical mechanics problem 
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in question would be identical to that of a ferromagnetic 
Ising model with 1/r 2 interactions in zero magnetic field. 
Similar to the Ising model there could be two phases: 

In the ordered (broken symmetry) phase the system, 
initially prepared in one of the minima of the potential 
V(6 C , 6 S ), will reside in there indefinitely long despite oc- 
casional excursions induced by thermal fluctuations into 
neighboring minima followed by returns into the origi- 
nal minimum. The excursion is a nucleus of neighbor- 
ing ground state constricted by two domain walls having 
opposite signs of the change A9 C . For the Ising model 
with 1/r 2 interactions the domain walls are known to 
be responsible for the destruction of order<22; the same is 
expected to hold for the model ([7]). Domain walls have 
their quantum-mechanical counterparts in instantons - 
extremal solutions of the equations of classical dynamics 
in imaginary time satisfying the condition that as imagi- 
nary time r traverses from minus infinity to infinity, the 
instanton interpolates between degenerate ground states. 

The similarity to the Ising model also implies that 
the broken symmetry phase of (7]) may be removed by 
proliferating domain walls. As a result the system will 
freely hop through the minima of V(9 C , 9 S ) leading to the 
destruction of the Coulomb staircase in the quantum- 
mechanical context. 

In order to understand quantitatively the role of the 
domain walls we find it useful to rewrite the action J7]) 
in an equivalent form which abandons the Fourier repre- 
sentation of the first term and shows explicitly that the 
latter vanishes for the uniform system: 

Aeff = -H 9* 1 f 0,(r)dTln(\T - T'\A )e\{T')dT' 

v=c,s J 

+ J drV{B c ,9 s ) (57) 

where the integration in the first term is limited by the 
condition \t — t'\Aq > 1 and the dots indicate derivatives. 

The domain wall of the field theory f57j) can be 
found by solving equations 5A e ff /S8 C}S — subject to 
the boundary conditions 6 C (— oo) — 9 c (oo) = ±q and 
9 S (— oo) — s (ao) = ±1. In the limit of large barriers 
the solution will have a form of a steep rounded step 
which can be approximated by 

4(r) = ±q5(r - t ), 6 s {t) = ±S(t-t ) (58) 

where To is arbitrary position of the domain wall. The 
combination of Eq . <f57|l and (|58|) implies that in the large- 
barrier limit the physics is dominated by thermally ex- 
cited logarithmically interacting domain walls whose se- 
quence is constrained by the zigzag locations of classical 
ground states discussed earlier. 

A. The "entropy" argument of Landau, Thouless 
and Kosterlitz and Thouless 

We begin with a semi-quantitative argument similar to 
the "entropy" arguments of Landau 2 ^, Thouless^, and 



Kosterlitz and Thouless^. The energy of a single do- 
main wall diverges with the system size L, as indicated 
by Eqs.JlZP and ([58]): E = h{q 2 g~ 1 + g' 1 ) m(A L). The 
entropy is given by the logarithm of the number of ways 
to put the domain wall on the line: S = ln(Ao-L). There- 
fore we find for the change of the free energy F = E — TS 
for the introduction of a single domain wall 

F = (h(^- + -)-t) ln(Aoi) (59) 

V 9c g s J 

The free energy changes sign at the temperature 

T c = + — ) (60) 

9c 9s 

This suggests that above T c the ground-state is unstable 
against spontaneous creation of domain walls. Below T c 
there are only bound pairs of domain walls while above 
T c there are free domain walls. 

To see the implications of these conclusions for the 
problem of a nearly-closed Coulomb-blockaded quantum 
dot we note that in the mapping between quantum 
field-theoretical models and those of classical statistical 
mechanics^ zero-point motion corresponds to thermal 
fluctuations, i. e. the role of Planck's constant h is played 
by the temperature T. Therefore provided 

1 q 2 

1 — < (61) 

9s 9c 

weak tunneling will not smear the Coulomb staircase. If 
this is the case, then there is non-zero lower limit to the 
magnitude of the population jump A9 C (given by q) at 
half-integer dimensionless gate voltage N: 

(A9 c ) rnm = (g c (l - g- 1 )) 1 ' 2 (62) 

This parallels the conclusion that there is non-zero lower 
limit to the magnetization of the one-dimensional clas- 
sical Ising model with 1/r 2 interactions^ 2 -. A conclusion 
qualitatively similar to l(62|) has been found earlier for 
the spinless version of the Coulomb blockade problem 7 . 
One may argue that a formal transition from spinfull to 
spinless electrons would consist in imposing SU(2) spin 
symmetry, g s = 2, followed by the replacement^ g c = 2g. 
Then Eq. ([62j) will exactly reproduce the spinless result. 

In the presence of ST/ (2) spin symmetry, g s = 2, the 
criterion (|6TT) becomes g c < 2q 2 which for nearly closed 
dot, q — 1, further simplifies to g c < 2. The latter con- 
dition is satisfied in the presence of short-range repul- 
sion between the electrons^. We conclude that repul- 
sive interactions between the electrons aid the Coulomb 
staircase to survive disordering effect of the zero-point 
motion. This meets the physical expectation as repul- 
sions decrease quantum fluctuations in the charge chan- 
nel, thus making the system more classical. 

It is curious that non-interacting electrons, g c>s = 2, 
and a nearly-closed dot, q = 1, correspond to the sign 
of equality in Eq. ([6Tj) - the "entropy" argument predicts 
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neutral stability. This is the first indication that the free- 
electron case has to be treated with special care. Here 
the analysis has to be conducted along two directions: 

First, renormalization-group analysis of Section IVB 
improves on the "entropy" argument, thus allowing us 
to make definite conclusions about the case of non- 
interacting electrons. Second, one should try to go be- 
yond the two-state approximation which is in the heart 
of the zigzag nature of the classical ground states. This 
more complicated problem will be left for future study. 

In the case of inequality opposite to lf6lj) . even weak 
tunneling destroys the Coulomb blockade. 

B. Renormalization-group treatment 



A 




One can improve on the "entropy" argument via a 
renormalization-group analysis of the action ([57| pertur- 
bative in dimensionless tunneling amplitude A (measured 
in units of M_o) between nearest ground states. Another 
pertinent parameter is the distance q between the near- 
est wells of the potential V(9 C ,9 S ) in the 9 C direction. 
Corresponding lowest order renormalization-group flow 
equations are implicit in the work of Kane and Fisher—: 



dA 
Hi 



(1 



1 



)A 



dl 



= -8q 2 A 2 



(64) 



If we set ■ g c — g s — ig, then Eqs.([63j and (J64J) reduce 
to those of Kane and Fished describing the physics of 
resonant tunneling of interacting spinless electrons onto 
and off a dot formed by two barriers in the presence of a 
gate which couples to the electrons between the barriers. 

On the other hand, setting g s = 2 and g c — 2g nearly 
reproduces renormalization-group equations describing 
spinless version of the problem of nearly-closed Coulomb- 
blockaded quantum dol£; the physics encoded in Eqs. l(63|) 
and ([64)1 parallels in many ways that of the spinless case. 

Comparing Eqs. ([6"Tj) and l|63p we conclude that the "en- 
tropy" argument of Section IVA is equivalent to first- 
order renormalization-group calculation establishing the 
range of stability of the zero-tunneling fixed line. 

Eqs. f63|) and |64|) have to be solved subject to the "ini- 
tial" conditions, A(l = 0) = Ao and q(l = 0) = go; the 
values of Ao and go can be computed from the form of 
the action (|57|l . 

The flow diagram corresponding to Eqs. f63|) and f64|) 
is sketched in Fig. 2. As typical of the Kosterlitz-Thouless 
phase transition 2 ^, the flows can be divided into two main 
regions - presently they correspond to the loci of the mod- 
els where the Coulomb blockade either survives disorder- 
ing effect of zero-point motion or ceases to exist. 

In the former case the systems satisfying 1 — gj 1 — 
Qo/gc < (compare with Eq.lJBTJ) and 4 A 2 , < q 2 /g c + 
fJs 1 - 1 - (1 - g^)bx(q 2 /g c {l - g' 1 )), i. e. those to 



Figure 2: The flow diagram corresponding to Eqs. ll63ll and 
|1C4[| . The stable part of the A = fixed line is shown in bold, 
and the arrows indicate the direction of the flow. The models 
whose loci are given by the "initial" conditions "a"-"e" all ex- 
hibit smeared Coulomb staircase with characteristic flow pa- 
rameter l c , given by Eqs. lf65l) - (l69l) , respectively. Specifically, 
non-interacting electrons, g c .s = 2, entering a nearly-closed 
dot, qo = 1, are represented by the "initial" condition 'd". 



(63) the right of the separatrix 4 A 2 = q 2 /g c + g s 1 — 1 



(1 — g s 1 ) ln(q 2 /g c (l — g s 1 )) are carried by the flow to 



the stable fixed line, A = 0, 1 — g s 



q 2 /g c < o. 



This indicates that tunneling is irrelevant in macro- 
scopic limit, and the discontinuity of the dot popula- 
tion, A8 C (given by g), at half- integer N satisfies the 
inequality {A9 c )min > (<?c(l — S^ 1 )) 1 ^ 2 - The equality is 
reached at the Kosterlitz-Thouless phase transition point, 
4A 2 C = q 2 /g c + g- 1 - 1 - (1 - g^ 1 ) ln( % 2 / 9c (l - g- 1 )), 
where the population jump ([62|l depends only on inter- 
particle interactions. This universal relationship which is 
a general feature of a Kosterlitz-Thouless transitional is 
analogous to that between the magnetization jump, the 
amplitude of the 1/r 2 interaction and the phase transi- 
tion temperature of the Ising model—,—. 

The systems to the left of the separatrix 4 A 2 = q 2 /g c + 
g^ 1 — 1 — (1— 57 1 ) m (<Z 2 / 5c (1— 97 1 )) are taken by the flow 
outside of the range of applicability of Eqs. f63|) and (|64| . 
the effective tunneling amplitude, A, is seen to increase 
while q, the distance between the minima of the potential 
V(9 C , 9 S ) in the 9 C direction, tends to zero. We are thus in 
the regime of smeared Coulomb staircase with continuous 
dependence of the dot population on the gate voltage. Al- 
though the Kosterlitz-Thouless transition is continuous, 
in going from sharp to smeared Coulomb staircase, the 
shape of the staircase changes dis continuously. 

In the classical limit of closed dot, Ao = 0, the dis- 
tance between the minima of the potential V(9 C ,9 S ) in 
the 9 C direction is unity, go = 1- Then SU(2) symmetric 
models with electron-electron repulsion, i. e. those sat- 
isfying g s — 2, g c < 2, belong to the stable part of the 
fixed line, A = 0, 1 — ,g s _1 — g 2 / g c < 0. Allowing weak 
tunneling, i. e. sufficiently small Aq, only renormalizes 
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the population jump at half-integer gate voltage N down- 
ward so that < A8 C < 1. This supports the conclusion 
of Section IVA that weak tunneling does not destroy the 
Coulomb blockade in the presence of interparticle repul- 
sions. However at sufficiently large tunneling amplitude 
A 0c the Kosterlitz-Thouless transition takes place, and 
for A > A 0c the Coulomb staircase ceases to exist. 

In this same classical limit, Ao = 0, qo = 1, but with 
the electrons in the channel satisfying 1— — <7 t T 1 > 

0, we find ourselves on the unstable part of the A = 
fixed line. Here allowing infinitesimally small tunneling, 
Ao, immediately destroys the classical staircase which is 
consistent with the "entropy" argument. 

The free-electron case, g c , s = 2, is marginal because 
classically, Ao = 0, qo = 1, it corresponds to the end 
point of the zero-tunneling, A = 0, stable fixed line. 
However allowing infinitesimally small tunneling, as in- 
dicated in Fig. 2, immediately moves it into the regime 
where the tunneling amplitude grows under renormaliza- 
tion, and thus the Coulomb blockade is destroyed. This 
is in agreement with the original result. 

The systems where the Coulomb blockade is destroyed, 

1. e. those which under renormalization transformation 
are taken outside the range of applicability of perturba- 
tive, A < 1, equations lf63|) and (f64|) . can be charac- 
terized by the characteristic flow parameter l c for which 
A(/ c ) ~ 1. The corresponding frequency scale, A c = 
Aoe _ic (see Eq. (|50j) ) marks the onset of strong charge 
(electron number) fluctuations on large time scales: 

At times shorter than A" 1 , the tunneling remains 
weak, A(Z) <C 1, and the dot population is frozen at an 
integer value. On the contrary, at times larger than A^ 1 , 
the tunneling is strong and the population of the dot fluc- 
tuates freely (around a half-integer value) thus making 
the particle discreteness irrelevant. The "spin" counter- 
part of the parameter l c (see Eq. (J52J) plays a related role 
in marking the onset of freezeout of spin fluctuations per- 
tinent to understanding the physics of nearly-open dot. 

The estimate of l c cannot be expected to be more than 
qualitative as the perturbative equations (|63| and f64|) 
are at most only indicative of its presence: 

If classically the model in question belongs to the un- 
stable part of the zero tunneling fixed line, then the flow 
starts out vertically away from the A = axis. If addi- 
tionally the beginning point is not too close to the end 
point A = 0, 1 — g^ 1 — q 2 g~ 1 — (initial condition "a" in 
Fig. 2), then to lowest order in Ao, the renormalization of 
q given by Eq. ljMf can be neglected. With this in mind 
the characteristic flow parameter l c can be found to be 



lr 



1 



In 



1 

Ao" 



(65) 



1 - 9s - g c 

If, on the other hand, the beginning point of the flow is 
close to the end point A = 0, 1 — g^ 1 — q 2 g c ~ 1 = (initial 
condition "b" in Fig. 2), then renormalization of q, given 
by Eq. (|64|) is no longer negligible; instead we find 



In 



9s 



9c 



(1 - 9 S 



Aoil-gT 1 ) 1 '* 



(66) 



We note that the estimates l(65j) and l(66j) are only loga- 
rithmically sensitive to the uncertainty, A(Z C ) ~ 1, built 
into the definition of the characteristic flow parameter l c . 

If the beginning point of the flow belongs to the out- 
going leftward separatrix in Fig. 2 (initial condition "c"), 
then to lowest order in A the characteristic flow param- 
eter l c is given by 



L 



1 



2A (2(l-. 9s 7 1 )) 1 / 2 



(67) 



while if the initial point is at the minimum of the flow 
trajectory, 1 — g^ 1 — q 2 g^ 1 = (initial condition "d" in 
Fig. 2), then we find 



lr 



4A (2(l-. 9s 7 1 )) 1 /2 



(68) 



The models located between the initial conditions "c" and 
"d" in Fig. 2 have characteristic flow parameter l c be- 
having as A " 1 (l — g^ 1 )^ 1 ^ 2 and interpolating between 
the estimates lf67j) and (j68"l) . We note that the systems 
represented by the initial condition "d" in Fig. 2 include 
non-interacting electrons, g c>s — 2, for which the charac- 
teristic flow parameter is l c = 7r/4Ao. 

If the model in question lies close to the point 
of Kosterlitz-Thouless phase transition separating the 
regimes of sharp and smeared Coulomb staircase, i. e. 
slightly higher than the incoming rightmost separatrix 
(initial condition "e" in Fig. 2), then the effective tun- 
neling rate A(l) initially decreases with I followed by an 
increase at larger I. As a result the characteristic flow 
parameter 



lr 



7T 



(69) 



4- v /(l- ffs - i )A 0c (A - A 0c )) 



is dominated by the "time" spent in the vicinity of the 
end point A = 0, 1 — g^ 1 — q 2 g^ x = 0. Here Ao c is 
the tunneling amplitude at the phase transition and it is 
assumed that Ao — Ao c -C Ao c . As expected, the corre- 
sponding scale in imaginary time direction L c = A ~ 1 e' c 
has a functional form characteristic of the behavior of 
the correlation length in the vicinity of the Kosterlitz- 
Thouless phase transitional. 

The curious feature of the leading order estimates 
(|67ll - l(69|) is that they are insensitive to the uncertainty 
A(Z C ) ~ 1 built into the definition of the characteristic 
flow parameter l c . 



C. Vicinity of half-integer dot population: 
heuristic argument 

One of the most intriguing statements made by 
Matveev is the argument that (for non-interacting elec- 
trons) the weak-tunneling results are implied by those 
found for a nearly-open dol£. 

Renormalization-group analysis of Section IVB enables 
us to see how and why the types of non-analytic behavior 
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of the ground-state properties found for a nearly-open dot 
(Section III) may be applicable well beyond that regime: 

This is expected to be the case for every system which 
under renormalization transformation flows outside the 
range of applicability of perturbative equations (|63|) and 
([64]) . i. e. into the region of large tunneling amplitude 
A and vanishing distance q between the minima of the 
potential V(9 C ,9 S ), Eq.(7])- Since this is the regime of a 
nearly-open dot, it seems plausible that no new physics 
intervenes for intermediate tunneling. As far the the 
weak-tunneling regime goes, the description of Section III 
must take over at frequencies smaller than A c = Aoe~' c . 

As a test of this statement we now show how the 
weak-tunneling results can be guessed from their strong- 
tunneling counterparts. First we remind the reader that 
as charge fluctuations renormalize the distance q between 
the zigzag minima of the effective potential V(9 Cl 9 s ) in 
the 9 C direction, the step of the zigzag in the 9 S direc- 
tion remains unaffected. The latter periodicity becomes 
"invisible" for q = and half-integer dimensionless gate 
voltage because the effective potential ceases to be de- 
pendent on the spin field 9 S . 

The periodicity with respect to the spin field 9 S be- 
comes "visible" again if the dimensionless gate voltage 
N is tuned away from half-integer value. For small de- 
viation the effective potential will have be proportional 
to SN, and the overall energy scale of the potential will 
be given by the charging energy Ec. These arguments 
imply the estimate 



V(0 S ) ~ E C 6N cos it0 s 



(70) 



where in displaying the 9 S periodicity we restricted our- 
selves to the leading Fourier harmonic and set zero of V 
at SN — 0. The potential of exactly the same functional 
form enters the effective action describing the strong- 
tunneling dynamics of the spin degrees of freedom (fT4l) , 
and the vicinity of a half-integer dimensionless gate volt- 
age, SN <C 1, guarantees the applicability of the results 
of Section III. Keeping in mind that the role of the fre- 
quency cutoff Ao will be played by A c = Aoe _ic , we find 
the correspondence 



E C SN 



-SNe 



(71) 



where we assumed a fixed ratio of the charging energy 
Ec to the original cutoff energy hA . Therefore the 
ground-state properties of nearly-closed quantum dot in 
the vicinity of half-integer dimensionless gate voltage can 
be extracted from the strong-tunneling results (Section 
III) provided the substitution (|7T|) is made and the fre- 
quency cutoff Ao is replaced with A c = Aoe~' c . 

For example, making these changes in Eq. ([29|) we find 
that in the presence of 577(2) spin symmetry, g s = 2, the 
non-analytic correction to the ground-state energy in the 
weak-tunneling regime has the form 



SE ~ ~E c e l "SN 2 In 



/ 1 



\SN 2 



(72) 



For non- interacting electrons, when l c = tt/AAq (see 
Eq. ([68|) ). and to leading order in the tunneling ampli- 
tude, this reproduces the functional form derived earlier 4 
via the analogy to the two-channel Kondo problem and 
the exact solution of the latter. 

Our analysis implies that non-analytic behavior exhib- 
ited in Eq. |72|) is also characteristic of any system with 
SU(2) spin symmetry having not very strong interparti- 
cle repulsions and sufficiently large tunneling amplitude. 
Specifically, in the vicinity of the critical tunneling am- 
plitude Ao c which corresponds to the Kosterlitz-Thouless 
transition into the regime of sharp Coulomb staircase, 
the correction to the ground-state energy is given by the 
combination of Eqs.((69]) and (|72"1) . 



V. CONCLUSIONS 

In this paper we analyzed the zero-temperature dy- 
namics of electrons hopping onto and off a large quan- 
tum dot connected to a reservoir via a one-dimensional 
channel and interparticle interactions inside the channel 
as well as Coulomb interactions among the electrons on 
the dot were taken into account. The main analytical 
tool applied to understand both the limits of strong and 
weak tunneling was the renormalization-group method. 

For strong tunneling (nearly-open dot) the charge 
discreteness is irrelevant and the physics is dominated 
by fluctuating spin degrees of freedom which are re- 
sponsible for non-analytic behavior of the ground-state 
properties at half-integer dot population. Specifically, 
the previously found logarithmic divergence of the dot 
capacitance^ is demonstrated to be a consequence of the 
assumption that electrons in the are non-interacting. We 
show that the same conclusion holds for physical systems 
possessing 577(2) spin symmetry. More generally, how- 
ever, the logarithmic singularity is replaced by a power 
law non-analyticity whose exponent is determined by the 
degree of spin fluctuations. 

For strong tunneling (nearly-closed dot) and half- 
integer dot population the spin degrees of freedom play 
a secondary role and the physics is dominated by zero- 
point hops of the electrons onto and off the dot. If these 
fluctuations are sufficiently strong, the Coulomb block- 
ade is destroyed by arbitrarily weak tunneling. Non- 
interacting electrons belong to this regime which is in 
agreement with previous studies^. If the charge fluctu- 
ations are sufficiently weak, then the Coulomb blockade 
survives provided the tunneling is not too strong. The 
regimes of sharp and smeared Coulomb staircases meet 
at a phase transition of the Kosterlitz-Thouless kind. 

Renormalization-group analysis of the weak-tunneling 
regime also makes a case in favor of the argument (first 
given for non-interacting electrons^) that non-analytic 
behavior of the ground-state properties at half-integer 
dot population found in the limit of nearly-open dot can 
extend well beyond the regime of strong tunneling. 

Finally, renormalization-group treatment of weak tun- 
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neling rests on the two-state approximation in which all 
the population states of the quantum dot separated from 
the two degenerate minima by energy gaps are ignored. 
Although physically plausible, similar to the spinless ver- 
sion of the problem- , its applicability to the case of non- 
interacting electrons requires further justification. In- 
deed, in the classical limit of closed dot, Ao = 0, qo = 1, 
the locus of non-interacting electrons, g c . s = 2, coin- 
cides with the end point of the zero-tunneling fixed line 
in Fig. 2. Allowing small finite tunneling A brings this 
model into the unstable regime characterized by the ini- 
tial condition "d". If we now assume that the ignored 
metastable states are populated with small finite proba- 
bility, then the actual electron number fluctuations will 
be suppressed stronger than predicted by the two-state 
approximation. 

Heuristically, the effect of weaker zero-point motion 
can still be seen within the two-state analysis - the ini- 
tial point of the flow will have to be located somewhat to 
the right of the initial condition "d" in Fig. 2. If the initial 



condition lies to the right of the incoming separatrix of 
the flow, then the Coulomb blockade can survive weak 
tunneling. Additionally, the Coulomb staircase would be 
destroyed at finite tunneling via a phase transition. This 
possibility calls for an investigation which avoids the two- 
state truncation altogether. Variational analysis of this 
problem which does not rely on the two-state approxi- 
mation and similar to the one used in the spinless case^ 
will be presented elsewhere^!. 
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